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This  paper  presents  an  adaptive  supervisory  control  strategy  for  a  fuel  cell/battery-powered  city  bus  to 
fulfill  the  complex  road  conditions  in  Beijing  bus  routes.  An  equivalent  consumption  minimization  strat¬ 
egy  (ECMS)  is  firstly  proposed  to  optimize  the  fuel  economy.  The  adaptive  supervisory  control  strategy 
is  exploited  based  on  this,  incorporating  an  estimating  algorithm  for  the  vehicle  accessorial  power,  an 
algorithm  for  the  battery  charge-sustaining  and  a  Recursive  Least  Squares  (RLS)  algorithm  for  fuel  cell 
performance  identification.  Finally,  an  adaptive  supervisory  controller  (ASC)  considering  the  fuel  con¬ 
sumption  minimization,  the  battery  charge-sustaining  and  the  fuel  cell  durability  has  been  implemented 
within  the  hybrid  city  buses.  Results  in  the  “China  city  bus  typical  cycle”  testing  and  the  demonstrational 
program  of  Beijing  bus  routes  are  presented,  demonstrating  that  this  approach  provides  an  improvement 
of  fuel  economy  along  with  robustness  and  ease  of  implementation.  However,  the  fuel  cell  system  does 
not  leave  much  room  for  the  optimal  strategy  to  promote  the  fuel  economy.  Benefits  may  also  result  in  a 
prolongation  of  the  fuel  cell  working  life,  which  needs  to  be  verified  in  future. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

With  growing  concerns  on  energy  crises  and  environmental 
issues,  Proton  Exchange  Membrane  (PEM)  fuel  cell  is  favored  for 
automotive  applications  because  of  zero  emission,  high  efficiency 
and  low  noise.  However,  it  has  disadvantages  of  high  price  and  short 
life-time,  which  makes  the  commercialization  difficult.  The  disad¬ 
vantages  could  be  partly  overcome  by  hybridization  [1].  There  are 
two  potential  benefits  by  combining  a  PEM  fuel  cell  with  an  energy 
storage  source.  Firstly,  the  durability  of  the  fuel  cell  stack  could 
be  improved  because  the  additional  energy  source  can  fulfill  the 
transient  power  demand  fluctuations.  Secondly,  the  ability  of  the 
energy  storage  source  to  recover  braking  energy  enhances  the  fuel 
economy  greatly  [2]. 

The  hybrid  configuration  raises  the  question  of  energy  man¬ 
agement  strategy,  to  achieve  system-level  performance  objectives 
such  as  fuel  economy,  battery  charge-sustaining,  fuel  cell  durability 
while  satisfying  the  system  constraints  [3]. 

The  recent  literature  contains  several  energy  management 
strategies  for  Hybrid  Electric  Vehicles  (HEV).  Energy  management 
strategies  have  been  investigated  for  both  internal  combustion 
engine  hybrid  vehicles  and  fuel  cell  hybrid  vehicles,  which  shows 
that  there  are  many  common  characteristics.  These  strategies  can 
be  categorized  into  three  types:  rule-based  strategies,  instanta¬ 
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neous  optimization  strategies  and  global  optimization  strategies. 
A  rule-based  strategy  can  be  easily  implemented  for  the  real-time 
applications  based  on  heuristics.  It  can  be  further  improved  by 
extracting  optimal  rules  from  optimal  algorithms  [4-6].  An  instan¬ 
taneous  optimization  strategy  evaluates  the  electrical  energy  from 
the  battery  using  an  equivalent  fuel  consumption  concept  [7-9]. 
A  global  optimization  strategy  minimizes  a  cost  function  over  a 
determined  driving  cycle  [10,11].  The  global  strategy  could  not  be 
implemented  in-vehicle,  but  it  can  be  regarded  as  a  basis  of  com¬ 
parison  for  the  evaluation  of  other  real-time  control  strategies. 

However,  still  there  are  some  special  characteristics  for  fuel  cell 
hybrid  vehicles.  The  fuel  cell  voltage  drops  with  increasing  current. 
A  reactant  starvation  occurs  at  high  currents  and  dynamic  loads 
because  the  transport  of  reactant  gases  is  not  able  to  keep  pace 
with  the  amount  used  in  the  reaction.  Thus,  the  fuel  cell  current 
and  the  dynamic  of  the  load  need  to  be  limited  [12].  Moreover, 
performance  degradation  of  the  PEM  fuel  cell  affects  the  battery 
charge-sustaining  and  the  vehicle  drivability.  An  online  identifi¬ 
cation  algorithm  is  therefore  necessary  to  identify  the  fuel  cell 
performance.  Control  parameters  should  be  adjusted  adaptively 
in  the  vehicle  supervisory  controller  to  achieve  an  overall  optimal 
performance. 

This  paper  deals  with  an  adaptive  supervisory  control  strat¬ 
egy  of  a  fuel  cell/battery-powered  hybrid  city  bus.  It  proceeds  as 
follows:  Section  2  describes  the  hybrid  powertrain  and  the  city 
bus.  Section  3  details  the  battery  equivalent  hydrogen  model  and 
the  equivalent  consumption  minimization  strategy  (ECMS).  The 
adaptive  supervisory  control  strategy  based  on  the  ECMS  is  pre¬ 
sented  in  Section  4,  incorporating  an  estimating  algorithm  for  the 
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Nomenclature 

a ,  b  fit  coefficients 

bfC  characteristic  coefficient  of  the  fuel  cell  polarization 

curve  (V) 

c  fit  coefficient 

C  hydrogen  consumption  of  the  hybrid  fuel  cell  pow¬ 

ertrain  (gs_1) 

Cbat  battery  equivalent  hydrogen  consumption  (g  s-1 ) 
Cfc  fuel  cell  hydrogen  consumption  (g  s-1 ) 

Qc.avg  average  fuel  cell  hydrogen  consumption  (g  s-1 ) 
d,f  fit  coefficients 

h  custom  vector 

/dc,tg  DC/DC  target  current  (A) 

Jfc  PEM  fuel  cell  net  output  current  (A) 
k  coefficient  of  the  energy  management  strategy  ( kW) 

I<  compensation  coefficient  for  battery  charge- 

sustaining 

K  coefficient  matrix 

K\  custom  coefficient 
P  coefficient  matrix 

Pbat  battery  power  (kW) 

^dc.avg  average  DC/DC  output  power  (kW) 

Pdc0  fit  coefficient  (kW) 

Pfc  fuel  cell  net  power  (kW) 

Rchg  battery  charging  resistance  (£2) 
fidis  battery  discharge  resistance  (Q) 
fifC  characteristic  coefficient  of  the  fuel  cell  polarization 

curve  (Q) 

Paux  vehicle  auxiliary  power  (kW) 

Pbat, opt  battery  optimal  power  (kW) 

Pm  electric  motor  actual  output  power  (kW) 

Pm,tg  electric  motor  target  power  (kW) 

Pmb.tg  electric  motor  target  brake  power  (kW) 

Pmd.tg  electric  motor  target  drive  power  (kW) 

Q.  battery  capacity  (A  h) 

SOC  state  of  charge 
SOCH  upper  limit  of  SOC 
SOCl  lower  limit  of  SOC 
SOCtg  battery  target  SOC 

SOCi  battery  SOC  up  to  the  energy  management  strategy 
SOC2  battery  SOC  up  to  the  vehicle  accessorial  power  esti¬ 
mating  algorithm 

SOC3  battery  SOC  up  to  the  target  drive  power  of  the  motor 
Pmd.tg 

SOC4  battery  SOC  up  to  the  target  brake  power  of  the 
motor  Pmb.tg 

Tqttg  electric  motor  target  torque  (Nm) 

Vbus  bus  voltage  (V) 

^bus.max  maximal  bus  voltage  (V) 

^bus.min  minimal  bus  voltage  (V) 

VfC  fuel  cell  terminal  voltage  (V) 

Vfc0  characteristic  coefficient  of  the  fuel  cell  polarization 

curve  (V) 

Vocv  battery  open  circuit  voltage  (V) 
xmin*  *max  custom  parameters 

Greek  symbols 

rjc hg.avg  average  charging  efficiency 
^dis.avg  average  discharging  efficiency 
rjm  motor  efficiency 
rdc  time  constant  for  the  low-pass  filter  (s) 
raux  time  constant  for  the  vehicle  accessorial  power  esti¬ 

mating  (s) 


/x,  k  custom  parameters 

ASOC  SOC  deviation  from  the  target 

0  custom  vector 

0  estimated  value  of  0 

Ardcl,  Ardc2,  Ardc3  modified  time  constants 


vehicle  accessorial  power,  an  algorithm  for  the  battery  charging- 
sustaining  and  a  Recursive  Least  Squares  (RLS)  algorithm  for  the 
fuel  cell  performance  identification.  Results  of  the  city  bus  in  the 
“China  city  bus  typical  cycle”  testing  and  in  the  demonstrational 
program  of  Beijing  are  described  in  Section  5.  Section  6  gives  the 
conclusions. 

2.  The  fuel  cell  hybrid  city  bus 

The  approach  presented  in  this  paper  was  implemented  in  sev¬ 
eral  city  buses,  which  were  used  as  the  “pick-up”  vehicles  in  the 
Marathon  Races  of  the  Beijing  Olympic  and  Paralympic  Games  in 
2008.  Fig.  1  shows  two  buses  on  the  Tiananmen  Square  on  24th 
August,  2008.  Each  of  them  is  equipped  with  ten  140L/20MPa 
hydrogen  tanks  and  weights  14,672  kg. 

As  electrochemical  devices,  fuel  cells  convert  chemical  energy 
into  electrical  energy  directly  without  mechanical  processes.  PEM 
fuel  cells  are  especially  preferred  in  automotive  applications 
because  they  are  efficient,  compact  and  have  a  good  ability  to  follow 
the  dynamic  loads.  The  output  power  of  a  single  cell  with  an  active 
area  of  200  cm2  is  less  than  100  W.  Therefore,  numerous  cells  are 
connected  in  series  to  form  a  multi-kW  stack,  and  a  system  contain¬ 
ing  several  stacks  generates  a  power  output  of  several  tens  of  kW. 
Besides,  a  fuel  cell  stack  also  contains  several  auxiliary  subsystems, 
such  as  air  supply  and  coolant  [2].  For  the  considered  city  buses, 
two  stacks  with  a  rated  power  of  40  kW  are  electrically  connected 
as  parallel  strings. 

A  Ni-MH  battery  has  the  advantage  of  good  charging/discharging 
performance  compared  to  a  Pb-acid  battery.  And  it  is  relatively 
cheap  compared  to  a  Li-ion  battery.  Thus,  a  Ni-MH  battery  with 
a  rated  capacity  of  80 Ah  is  utilized  for  the  considered  city  bus. 
The  battery  comprises  four  modules,  each  of  which  contains  70 
cells.  The  Ni-MH  cells  installed  in  the  vehicle  have  a  rated  volt¬ 
age  of  1.2  V.  Altogether,  the  280  cells  yield  a  rated  voltage  of  338  V. 
A  management  system  is  designed  and  mounted  on  the  battery 
to  equilibrate  the  cell  voltage,  avoid  the  overcharging  of  specific 
cells  and  measure  the  working  current,  the  terminal  voltage  and  the 
temperature. 

The  city  bus  is  powered  by  an  AC  motor  with  a  rated  power  of 
100  kW,  a  peak  power  of  150  kW  and  a  maximum  torque  of  1121  Nm. 
The  fuel  cell  stack,  the  Ni-MH  battery  and  the  AC  motor  are  con- 


Fig.  1.  Two  PEM  fuel  cell/battery  hybrid  city  buses  on  the  Tiananmen  Square. 
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Fig.  3.  The  time  triggered  distributed  control  system  based  on  TTCAN. 


nected  as  shown  in  Fig.  2.  Because  the  dc  bus  voltage  Vbus  is  lower 
than  the  fuel  cell  terminal  voltage  Vfc,  the  fuel  cell  stack  is  linked 
to  the  dc  bus  by  a  current-regulated  buck  converter,  which  delivers 
a  directional  current.  A  direct  current  to  direct  current  converter 
(DC/DC)  with  a  rated  power  of  80  kW  is  selected  to  mach  with  the 
PEM  fuel  cell  system. 

The  hybrid  vehicle  utilizes  the  PEM  fuel  cell  as  the  main  power 
source  and  the  battery  as  the  auxiliary  power  source.  The  current- 
regulated  DC/DC  is  used  to  control  the  output  power  of  the  fuel 
cell.  The  work  of  the  hybrid  powertrain  is  coordinated  by  a  super¬ 
visory  controller  within  a  TTCAN  (Time  Triggered  Controller  Area 
Network)  system,  which  is  shown  in  Fig.  3. 

The  TTCAN  network  communication  is  based  on  a  periodic  time- 
triggered  mechanism.  It  combines  the  advantages  of  event  and  time 
triggered  communications  in  distributed  real-time  control  systems. 
Main  advantages  of  time  triggered  networks  are  the  quasi  deter¬ 
ministic  behavior  and  ease  of  implementation  of  fault  tolerant 
systems.  One  of  the  disadvantages  is  the  restrictive  design  pro¬ 
cess,  because  the  communication  structure  needs  to  be  defined  in 
advance  [12,13]. 

The  time  triggered  distributed  control  system  comprises  a  vehi¬ 
cle  supervisory  controller  (VSC),  a  battery  management  system 
(BMS),  a  DC/DC  controller  (DCC),  a  motor  controller  unit  (MCU) 
and  a  fuel  cell  control  system  (FCC),  as  shown  in  Fig.  3.  Param¬ 
eters  for  the  vehicle  and  the  hybrid  powertrain  are  presented  in 
Table  1. 


Table  1 

Parameters  for  the  vehicle  and  the  hybrid  powertrain. 


Parameter  (unit) 

Value 

Vehicle  mass,  m  (kg) 

1.4672  x  104 

Frontal  area,  A  (m2) 

7.5 

Drag  coefficient,  CD 

0.7 

Rolling  resistance  coefficient 

o 

X 

00 

Mechanical  efficiency,  rjj  (%) 

95 

Mass  factor 

1.1 

PEM  fuel  cell  rated  power  (kW) 

80 

DC/DC  rated  power  (kW) 

80 

Ni-MH  battery  rated  capacity  (A  h) 

80 

Electric  motor  peak  power  (kW) 

150 

Electric  motor  peak  torque  (Nm) 

1.121  x  103 

Electric  motor  rated  power  (kW) 

100 

Electric  motor  maximal  rotational  speed  (rrnim1) 

6x  103 

3.  The  equivalent  consumption  minimization  strategy 
(ECMS) 


The  concept  of  equivalent  fuel  consumption  was  proposed  by 
Paganelli  et  al.  for  the  development  of  an  instantaneous  optimiza¬ 
tion  energy  management  strategy.  In  fuel  cell  hybrid  vehicles,  both 
hydrogen  energy  and  electrical  energy  can  be  used.  If  some  energy 
is  drawn  from  the  battery  at  the  current  sample  time,  the  battery 
will  have  to  be  recharged  to  maintain  the  state  of  charge.  The  energy 
will  be  provided  by  the  fuel  cell  system.  That  will  imply  extra  hydro¬ 
gen  consumption.  If  the  battery  is  recharged  at  the  current  sample 
time,  a  discharge  of  the  battery  is  required  to  maintain  the  state  of 
charge.  This  discharge  will  lead  to  a  reduction  of  the  hydrogen  con¬ 
sumption  in  the  future.  The  electrical  energy  consumption  of  the 
battery  is  transformed  into  an  equivalent  hydrogen  consumption 
to  make  the  two  comparable  [9,14]. 

The  output  power  of  the  fuel  cell  system  is  controlled  by  the 
current-regulated  DC/DC.  An  optimal  DC/DC  power  Pdcopt  is  cal¬ 
culated  to  minimize  the  hydrogen  consumption  of  the  hybrid 
powertrain  C,  which  comprises  the  actual  fuel  consumption  Cfc  and 
the  battery  equivalent  fuel  consumption  Cbat.  Aiming  at  minimiz¬ 
ing  the  vehicle  hydrogen  consumption,  the  equivalent  consumption 
optimization  problem  is  mathematically  formulated  as  follows: 


Pdc, opt  =  argminC  =  argmin(Cfc  +  /cCbat) 

Ale  Ajc 


(1) 


where  Pdc  is  the  DC/DC  output  power  (kW).  k  is  a  penalty  coefficient, 
which  biases  the  battery  equivalent  fuel  consumption  up  or  down 
depending  on  the  battery  SOC  (state  of  charge)  deviation  from  the 
target  [3]. 

The  optimization  problem  is  subject  to  several  constrains,  such 
as  battery  charging-sustaining  and  respecting  the  components’  lim¬ 
itations: 


subject  to  : 


'  SOCL  <  SOC  <  SOCH 

<  Vbus  ,min  <  14>us  <  Vbus  ,max 
^  0  <  Pdc  <  Pdc,max 


(2) 


where  SOCL  is  the  lower  limit  of  SOC,  SOCH  is  the  upper  limit  of  SOC, 
Vbus  is  the  dc-link  bus  voltage,  Vbus  max  and  Vbus  min  are  the  maximal 
and  minimal  allowed  value  of  Vbus  and  Pdc,max  is  the  maximal  output 
power  of  the  DC/DC. 

According  to  the  experimental  data,  the  fuel  cell  hydrogen  con¬ 
sumption  Cfc  can  be  written  as: 


f  t*Pdc  +  h  Pdc  >  Pdc0 

\  cPdc  +  d^dc  +  /  Pdc  <  PdcO 


where  a,  b,  c,  d,f  and  Pdc0  are  fit  coefficients. 

In  order  to  calculate  the  battery  equivalent  hydrogen  con¬ 
sumption  Cbat,  the  average  values  (average  charging/discharging 
efficiencies,  average  DC/DC  power  and  average  fuel  cell  hydrogen 
consumption)  are  used  because  the  operating  points  of  the  PEM 
fuel  cell  stack  and  the  battery  in  the  future  are  not  known  [9]. 
This  concept  is  further  developed  by  Fie  et  al.,  where  the  battery 
charging/discharging  efficiencies  are  calculated  based  on  the  Rint 
model  [15].  Then,  the  battery  equivalent  fuel  consumption  Cbat  can 
be  written  as: 


Q>at  — 


Pbat  Ffc.avg 
^chg.avg  Alc.avg 


Fbat  >  0 


Fbat^dis.avg 


Qc.avg 

Alc.avg 


4KchgPbat  \ 

y  ~^r) 


Fbat  <  0 


(4) 
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A  linear  penalty  coefficient  k  is  used  to  accomplish  the  charge- 
sustaining  operation.  It  is  defined  as: 


K=\  2/i 


(SOC  -  0.5(SOCh  +  SOCL)) 
SOCH  -  SOCL 


(5) 


where  /x  is  a  constant,  which  can  be  adjusted  to  reflect  the  battery 
charge  and  discharge  characteristics  [3]. 

For  the  hybrid  powertrain  under  discussion,  the  fuel  cell  +  DC/DC 
system  works  in  the  linear  area  during  the  most  of  the  time  on  the 
road.  Therefore,  the  optimized  problem  defined  in  Eqs.  (1)  and  (2) 
could  be  simplified  and  the  analytic  solution  to  the  problem  is  as 
follows: 


vbus  min(^ocv  ^bus.min) 
^dis 


A>at,opt  —  < 


,  K i  <  xm 


4Kdi 
0,  1  <  /Ci  < 


?  *min  <  ^1  <  1 

1 

^chg.avg^dis.avg 


^ocv(  ^  (^1  ^chg.avg^dis  ,avg) 

4^dis 

Vbus,max(Vbus,max  -  Vocv) 


Aic, opt  —  max 


^chg 

,tg(Wm 


^chg.avg^dis.avg 

Ki  >  ■  *m 


<  Ki  < 


^chg,avg^dis,avg 


^chg,avg^dis,avg 


(  ■  f  iTq  tgd 

min  — -i-5- 
\  \  Om 


+  Paux  ~  Anat  ,opt>  Pdc,r 


„),) 


(6) 


where  oom  is  the  electric  motor  rotational  speed  (rad  s-1 ),  ijm  is  the 
electric  motor  efficiency,  Paux  is  the  vehicle  auxiliary  power  (kW), 
Tqftg  is  the  electric  motor  target  torque  (Nm),  Pbatopt  is  the  battery 
optimal  power  (kW),  V0cv  is  the  open  circuit  voltage  of  the  battery 
(V)  and  K\ ,  xmin,  xmax  are  custom  defined  variables. 

The  plus  sign  “+”  is  used  when  the  electric  motor  works  as  a 
motor,  the  minus  sign  is  used  when  it  works  as  a  generator.  K\, 
xmin*  *max  are  defined  as: 


1  —  W  1  +  ^bus,i 


Vbu 


,-v0( 


V£* 


< 

*max 


'l+4Vbus,n 


^4>us,max  ^ocv 


V&v 


(7) 


K 1 


K 


^chg.avg 


^bus.min  and  l/bus,max  are  the  minimal  and  the  maximal  allowed 
bus  voltages  in  the  hybrid  vehicle.  V0cv  is  the  open  circuit  voltage 
of  the  battery,  which  is  a  function  of  the  battery  SOC.  Therefore, 
Pb at, opt  is  a  function  of  SOC  and  constant  /x. 

The  constant  /x  is  selected  to  balance  the  battery  SOC  over  a 
whole  trip.  For  the  considered  bus  and  the  routes  in  Beijing,  0.6  is  a 
suitable  value.  Fig.  4  indicates  the  relationship  between  Pbat.opt  and 
SOC  when  /x  =  0.6.  SOCtg  is  the  target  value  of  the  battery  SOC.  It  is 
defined  as  the  SOC  of  the  intersection  of  the  battery  optimal  power 
curve  and  the  x-axis. 

The  output  current  of  the  DC/DC  is  regulated  to  control  the  out¬ 
put  power  of  the  fuel  cell  stack.  A  low-pass  filter  is  used  to  calculate 
the  DC/DC  target  current  /dctg  in  order  to  avoid  the  fuel  cell  starva¬ 
tion  problem  under  dynamic  loads  as  follows: 


^dc.tg  — 


F’dc.opt/^b  us 
rdcs  + 1 


(8) 


where  rdc  is  the  time  constant  of  the  low-pass  filter  and  Pdc.opt  is 
calculated  as  in  Eq.  (6). 

The  ECMS  referred  above  is  the  core  of  the  adaptive  supervisory 
control  strategy.  Flowever,  it  is  not  enough  when  the  bus  is  driven 
in  bus  routes. 


Fig.  4.  Relationship  between  Pbat,opt  and  SOC  when  /x  =  0.6.  SOCtg  is  the  target  value 
of  the  battery  SOC.  It  is  defined  as  the  SOC  of  the  intersection  of  the  battery  optimal 
power  curve  and  the  x-axis. 


A  lag  between  the  applied  load  on  the  fuel  cell  and  the  response 
of  the  reactant  supply  system  results  in  an  undersupply  of  reactants 
to  the  fuel  cell,  which  leads  to  a  breakdown  of  the  chemical  reac¬ 
tion  and  a  rapid  loss  in  voltage  [2].  The  dynamics  of  the  load  need 
to  be  restricted  to  avoid  this  phenomenon.  In  the  ECMS,  the  load 
dynamics  are  restricted  by  the  low-pass  filter  with  a  time  constant 
rdc.  Flowever,  the  performance  of  the  fuel  cell  degrades  after  a  long 
working  time.  The  time  constant  should  be  changed  according  to 
the  status  of  the  fuel  cell  system.  An  online  identification  algorithm 
for  the  fuel  cell  performance  is  essential. 

What  is  more,  the  vehicle  auxiliary  power  faux,  which  is  con¬ 
sumed  by  the  accessorial  components,  e.g.  the  air  condition  and  the 
electric  power  steering  system,  changes  on  the  road.  The  supervi¬ 
sory  controller  should  have  the  ability  to  identify  the  power  online. 
Besides,  the  battery  may  be  over  charged  by  braking  energy  when 
there  are  too  many  red  lights  in  the  bus  route.  An  algorithm  for  the 
battery  charge-sustaining  is  also  necessary. 


4.  The  adaptive  supervisory  control  strategy  (ASC) 

The  adaptive  supervisory  control  strategy  is  developed  to  ful¬ 
fill  the  complex  road  conditions  in  bus  routes.  Based  on  the  ECMS, 
the  adaptive  strategy  includes  three  algorithms,  an  estimating  algo¬ 
rithm  for  the  vehicle  accessorial  power,  an  algorithm  for  the  battery 
charge-sustaining  and  a  Recursive  Least  Square  (RLS)  algorithm  for 
the  fuel  cell  performance  identification. 


4.1.  The  vehicle  accessorial  power  estimating  algorithm 


The  vehicle  accessorial  power  Paux  can  be  calculated  as: 


Paux  —  ^busUdc  +  4>at  —  An) 


(9) 


where  /dc  is  the  actual  output  current  of  the  DC/DC  (A).  It  is  mea¬ 
sured  by  the  DCC;  /bat  is  the  actual  output  current  of  the  battery 
(A).  It  is  measured  by  the  BMS;  Jm  is  the  actual  input  current  of  the 
motor  (A).  It  is  measured  by  the  MCU. 

Flowever,  the  result  of  Eq.  (9)  includes  extra  noise,  which  may  be 
caused  by  the  electromagnetic  interference  of  the  electric  compo¬ 
nents  in  the  hybrid  powertrain.  A  low-pass  filter  is  used  to  extract 
the  actual  value.  The  filtered  accessorial  power  could  be  expressed 
as: 


Faux  = 


^busUdc  +  fbat  ^m) 
TauxS  +  1 


(10) 
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Fig.  5.  The  original  signal  of  the  vehicle  accessorial  power  Paux  with  extra  noise  and 
the  filtered  Paux  with  a  low-pass  filter,  raux  =  100  s. 


where  raUx  is  the  time  constant.  Fig.  5  gives  an  example  of  the  fil¬ 
tered  accessorial  power  and  the  original  signal  with  extra  noise, 
raux  =  100s. 


4.2.  The  algorithm  for  the  battery  charge-sustaining 


The  battery  SOC  may  deviate  from  the  target  on  an  actual 
bus  route,  e.g.  a  bus  route  with  many  red  lights.  An  algorithm 
is  therefore  necessary  to  compensate  the  SOC  deviation  from  the 
target. 

Based  on  the  ECMS  in  Section  3  and  the  estimating  algorithm  in 
Section  4.1,  the  battery  SOC  could  be  sketched  as  in  Fig.  6.  It  can  be 
formulated  as: 


SOC  =  SOQ  +  SOC2  +  SOC3  +  SOC4  ( 11 ) 

SOCi  is  up  to  the  energy  management  strategy.  It  can  be  written 


as: 


SOCt  = 


fcSOC 


tg 


^ocvQTdc^  +  i/ocvQs  +  k 


(12) 


where  Q.  is  the  capacity  of  the  Ni-MH  battery,  k  is  a  coefficient  of 
the  ECMS.  It  is  defined  as: 


k  = 


Pb, 


at, opt 


SOC  -  SOCtg 

d  Pbat.opt 


SOC  +  SOC 


tg 


(13) 


dSOC 


,  SOC  =  SOC 


tg 


SOC2  is  up  to  the  vehicle  accessorial  power  estimating  algorithm. 
It  can  be  expressed  as: 


SOC2  = 


— F>aux(rcjcTaUxS2  +  (ldc  +  raux)s) 
(i/ocvQ't dcs2  +  1/ocvQs  +  k)(rauxs  +  1 ) 


(14) 


SOC3  depends  on  the  drive  power  of  the  motor  Pmdtg.  It  can  be 
calculated  as: 


SOC3  = 


_ Cic  /  hmS _ p 

VqcvQT dc^  +  VqcvQs  +  k  ’  £ 


(15) 


where  pm  is  the  motor  efficiency.  SOC4  is  up  to  the  brake  power  of 
the  motor  Pmb>tg.  It  can  be  formulated  as: 


SOC4  = 


—  TdchmS  —  hm 
VocvQX dc^2  +  V0CvQs  +  1< 


^mb.tg 


(16) 


The  steady  state  value  of  SOC/,  i  =  1  —4  when  t^+00  could  be 
calculated  as: 


k  SOC 


tg 


SOCi  lt-00  VocvQtdcS2  +  VocvQs  +  k 
SOC2l^oc  =  0 

SOC3  lt-^-00  =  0 

enr  i  “^mPmb.tg 

0UL4  |t^oo  =  - - - 


I  =  soc 


tg 


(17) 


Thus,  the  steady  deviation  of  the  battery  SOC  from  the  target 
SOCtg  is 


ASOC  =  SOC|t 


)  —  S0Ct2  — 


tg 


(18) 


In  order  to  compensate  this  deviation,  the  target  current  of 
DC/DC  should  be  modified  as: 


J^mPmb.tg  nqs 

Me, tg, modified  —  Mc,tg  +  y  _|_  /^  v 

where  I<  is  a  compensating  coefficient  and  /dc  tg  is  calculated  as  in 
Eq.(8). 

Theoretically  speaking,  ASOC  ->  0  when  t -+  +00  if  I<  =  1.  To  elim¬ 
inate  the  deviation  in  a  limited  time,  I<  is  always  larger  than  1, 
depending  on  the  actual  bus  route.  Fig.  7  compares  the  simulation 
results  of  the  battery  charge-sustaining  algorithm  with  I<=  0  (the 
fine  line),  I<=  1  (the  dashed  line),  I<=  2  (the  bold  line)  in  a  certain 
bus  route. 


4.3.  The  RLS  algorithm  for  the  fuel  cell  performance  identification 

The  output  power  of  the  fuel  cell  system  declines  after  a  long 
driving  distance.  The  performance  degradation  can  be  clearly  iden¬ 
tified  by  comparing  the  polarization  curves  at  the  beginning  and  at 


Fig.  7.  Comparison  of  the  simulation  results  of  the  battery  charge-sustaining  algo¬ 
rithm  with  different  I<,  I<=  0  (the  fine  line),  I<=  1  (the  dashed  line)  and  I<=  2  (the  bold 
line). 
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Fig.  8.  Relationships  between  the  modified  time  constants  and  the  parameters  for 
the  polarization  curve. 

the  end  of  a  driving  distance  [16].  The  relationship  between  the  fuel 
cell  terminal  voltage  VfC  and  its  output  current  JfC  can  be  noted  as: 

Vfc  =  h'0 

0  =  (VfcO,Rfc,bfc)J  (20) 

h  =  (l-/fc-ln(/fc))T 

where  h  and  0  are  vectors,  Vfc0,  RfC  and  bfc  are  the  characteristic  coef¬ 
ficients  of  the  polarization  curve  [17].  According  to  the  RLS  theory, 
the  vector  0  can  be  estimated  as: 

"  0(fe)  =  0(k-  1)4-  K(k)[  Vfc(k)  -  h(k)T0(k  -  1 )] 

<  K(k)  =  P(k-  l)h(k)[/i(fc)TP(k-  l)fc(fe)  +  1]_1  (21) 

.  P(.k)  =  [/  -  K(k)h(k)J]P(k  -  1 ) 

where  0  is  the  estimated  value  of  0 ,  K  and  P  are  coefficient  matrixes 
[  18  ].  The  parameters  for  the  fuel  cell  polarization  curve  at  the  begin¬ 
ning  of  the  long  driving  distance  on  the  testbench  are  used  as  the 


200  400  600  800  1000  1200  1400  1600 


t(  S) 


Fig.  9.  Simulation  results  of  the  RLS  algorithm:  (a)  the  three  characteristic  parame¬ 
ters  for  the  polarization  curve  and  (b)  modified  time  constant  rdc>modified- 


Driver 

request 


Fig.  10.  Implementation  of  the  adaptive  supervisory  controller  (ASC)  in  the  hybrid  city  bus. 
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Fig.  11.  Energy  flow  diagrams  in  the  “China  city  bus  typical  cycle”:  (a)  a  rule-based  strategy,  (b)  the  ACS  without  braking  energy  regeneration  and  (c)  the  ACS  with  braking 
energy  regeneration. 
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initial  values  for  the  algorithm.  They  are 

0(0)  =  (  542.5059  0.2337  21.1059  )T  (22) 


The  time  constant  rdc  is  then  modified  according  to  the  results 
of  the  RLS  algorithm  as: 


< 


^clc,  modified 

Atdcl  =  /l(Vfco) 
ATdc2  =/2(^fc) 
l  Ardc3  =/3(bfc) 


.  Ardcl  Ardc2  ArdC3 

rdc  H - 5 - 1 - 5 - 1 - 5 — 


(23) 


Ardcl,  Ardc2  and  Ardc3  are  the  modified  time  constants,  which 
are  functions  of  Vfc0,  RfC  and  bfc  respectively.  On  one  hand,  perfor¬ 
mance  degradation  of  the  fuel  cell  stack  leads  to  a  decrease  in  Vfc0, 
an  increase  in  Rfc  and  bfc.  On  the  other  hand,  a  larger  DC/DC  time 
constant  rdc  is  needed  to  make  the  output  power  of  the  fuel  cell 
stack  change  much  more  slowly  when  the  output  ability  declines. 
Therefore,  Ardcl  decreases  with  the  Vfc0  increase,  Ardc2  increases 
with  the  Rfc  increase,  Ardc3  increases  with  the  bfc  increase.  Fig.  8 
gives  an  example  of  the  three  curves. 

Fig.  9  gives  an  example  of  the  RLS  algorithm.  The  estimated 
parameter  converges  to  a  small  neighborhood  of  the  actual  value 
a  certain  period  after  the  system  initialization.  The  parameter  is 
set  to  the  initial  value  as  in  Eq.  (22)  when  the  estimated  result 
does  not  converge  to  the  neighborhood.  After  that,  the  estimated 
parameters  are  used  in  the  energy  management  strategy.  In  the 
given  example,  the  three  estimated  parameters  of  the  polarization 
curve  are  536  V,  0.38  £2  and  21.2  V  respectively.  The  convergence 
value  appears  600  s  after  the  system  initialization.  As  a  result, 
the  DC/DC  time  constant  is  adjusted  from  5  s  to  25  s,  which 
makes  the  dynamics  of  the  load  on  the  fuel  cell  much  smaller 
than  at  the  beginning.  Thus,  the  fuel  cell  stack  can  work  at  a 
quasi-steady  state.  And  the  working  life  of  the  fuel  cell  may  be 
prolonged. 


5.  Results  in  the  cycle  testing  and  the  demonstrational 
program  of  Beijing  bus  routes 

The  implementation  of  the  adaptive  supervisory  controller 
(ASC)  within  the  city  bus  is  described  in  the  control  scheme  shown 
in  Fig.  10.  It  is  a  synthetic  strategy  considering  the  optimal  system 
efficiency,  the  battery  charge-sustaining  and  the  fuel  cell  durabil¬ 
ity.  Several  city  buses  with  the  proposed  ASC  have  been  successfully 
demonstrated  in  the  marathon  races  of  the  29th  Olympic  and  13th 
Paralympic  Games.  Furthermore,  they  are  demonstrated  in  the  Bei¬ 
jing  bus  routes. 

The  fuel  cell  hybrid  city  buses  have  been  tested  in  the  “China  city 
bus  typical  cycle”  [19,20].  The  battery  SOC  was  maintained  around 
70%.  Fig.  ll(a)-(c)  shows  the  energy  flow  diagrams  with  a  rule- 
based  strategy,  with  the  ASC  and  with  the  ASC  +  braking  energy 
regeneration  respectively. 

The  energy  flow  diagram  with  a  rule-based  strategy  is  shown  in 
Fig.  11(a).  The  hydrogen  chemical  energy  is  calculated  on  the  basis  of 
the  LFIV  (Low  Fleat  Value).  Totally  67.313  MJ  (100%)  was  consumed 
during  the  cycle  testing.  Defined  as  (net  power)/(fuel  power),  the 
average  net  efficiency  of  the  fuel  cell  was  about  48.3%,  51.7%  of  the 
fuel  energy  was  lost  in  the  stack  and  its  auxiliary  components.  About 
7.14%  of  the  fuel  energy  was  consumed  by  the  vehicle  accessorial 
components.  The  average  efficiencies  of  the  electric  motor  were 
85%.  The  average  powertrain  efficiency  was  about  31.7%  and  the 
fuel  economy  was  9.5  kg  (100  km)-1. 

Fig.  11(b)  shows  the  energy  flow  diagram  with  the  ASC.  Totally 
61.383  MJ  (100%)  was  consumed  during  the  cycle.  The  average  net 
efficiency  of  the  fuel  cell  was  about  49.7%.  The  average  power¬ 


Fig.  12.  Fuel  cell  net  efficiency  (defined  as  (net  power)/(fuel  power))  vs.  fuel  cell  net 
power.  The  fitting  curve  and  the  experiment  data  for  the  considered  system. 


train  efficiency  was  about  34.9%  and  the  fuel  economy  was  9.3  kg 
(100  km)-1. 

Fig.  11(c)  shows  the  energy  flow  diagram  with  the  ASC  and  the 
braking  energy  regeneration.  Totally  50.072  MJ  (100%)  was  con¬ 
sumed.  The  average  net  efficiency  of  the  fuel  cell  was  about  49.5%. 
The  average  powertrain  efficiency  was  about  43.2%  and  the  fuel 
economy  was  7.8  kg  (100  km)-1 . 

Compared  to  the  rule-based  strategy,  the  ASC  and  the  braking 
energy  regeneration  lower  the  fuel  consumption  by  2%  and  16% 
respectively.  The  braking  energy  regeneration  contributes  much 
more  than  the  ASC  to  lowering  the  fuel  consumption.  This  is  because 
the  fuel  cell  system  does  not  leave  much  room  for  the  optimal  strat¬ 
egy  to  promote  the  fuel  economy.  Fig.  12  illustrates  the  relationship 
between  the  fuel  cell  efficiency  and  the  net  power.  The  fuel  cell  net 
efficiency  changes  slightly  between  20  kW  and  70  kW.  No  matter 
with  the  ASC  or  a  rule-based  strategy,  the  PEM  fuel  cell  follows 
the  vehicle  load  demand.  Hence,  there  is  not  much  room  left  for 
optimization. 

Fig.  13  presents  the  performance  degradation  of  the  fuel  cell  sys¬ 
tem  during  the  8000  km  mileage  in  the  demonstrational  program. 
Pfc  iooA  is  defined  as  the  net  output  power  of  the  fuel  cell  system 
at  a  net  output  current  of  100  A.  It  is  regarded  as  the  output  ability 


Fig.  13.  Performance  degradation  of  the  fuel  cell  system  during  8000  km  mileage. 
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of  the  fuel  cell  system.  The  fuel  cell  performance  declines  at  a  rate 
of  -0.26  W  km-1  (-0.628%  (1000  km)-1)  according  to  the  on-road 
data  in  the  8000  km  driving  distance.  Besides,  the  battery  SOC  was 
always  kept  around  70%  in  the  bus  routes. 

6.  Conclusions 

An  adaptive  supervisory  controller  has  been  exploited  consid¬ 
ering  the  fuel  economy,  the  battery  charge-sustaining  and  the 
fuel  cell  durability.  The  control  strategy  comprises  four  parts,  an 
ECMS,  a  vehicle  accessorial  power  estimating  algorithm,  a  battery 
charge-sustaining  algorithm  and  a  RLS  algorithm  for  the  fuel  cell 
performance  identification. 

The  analytical  solution  of  the  ECMS  is  derived  based  on  the  bat¬ 
tery  Rint  model.  The  battery  optimal  power  is  a  function  of  the 
battery  SOC  and  a  custom  defined  variable  /x,  which  can  be  adjusted 
to  reflect  the  battery  charge  and  discharge  characteristics.  The  SOC 
deviation  is  compensated  by  modifying  the  DC/DC  target  current 
based  on  the  motor  braking  power  and  a  compensating  coefficient 
I<.  The  variables  fi  and  I<  should  be  selected  according  to  the  actual 
bus  route. 

Compared  to  a  rule-based  strategy,  the  fuel  consumption  of  the 
city  bus  with  the  ASC  in  the  “China  city  bus  typical  cycle”  was  low¬ 
ered  from  9.5  kg  (100  km)-1  to  9.3  kg  (100  km)-1 .  It  can  be  lowered 
to  7.8  kg  (100  km)-1  with  the  ASC  +  braking  energy  regeneration 
strategy.  The  braking  energy  regeneration  contributes  much  more 
than  the  ASC  to  lowering  the  fuel  consumption.  This  is  because  the 
fuel  cell  system  does  not  leave  much  room  for  the  optimal  strategy 
to  promote  the  fuel  economy. 

During  the  8000  km  mileage  in  the  Beijing  bus  route  demonstra- 
tional  program,  the  output  power  of  the  fuel  cell  system  declines 
at  a  rate  of -0.26  W  km-1  (-0.628%  (1000  km)-1).  The  battery  SOC 
was  always  maintained  around  70%  in  the  bus  routes. 

The  proposed  approach  provides  an  improvement  in  fuel  econ¬ 
omy  along  with  robustness  and  ease  of  implementation.  Benefits 


may  also  result  in  a  prolongation  of  the  fuel  cell  working  life,  which 
needs  to  be  verified  in  future. 
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